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Abstract 


This article presents a new control chart for monitoring reliability using sudden death testing under the neutrosophic statistics 
(NS). The average run lengths of the in-control and the out-of-control process have been determined for evaluating the quick 
detection ability for small and moderate shifts. For the industrial use, tables and figures have been presented for different 
parameters. The proposed control chart is efficient in comparison with the existing control chart under classical statistics 
and value addition in the toolkit of the quality control personnel. 
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Introduction 


A control chart is an important tool of statistical process 
control for maintaining and improving the quality of the 
product to compete in the market [1]. An efficient control 
chart is one, which quickly indicates any special cause 
of variation in the process so that the engineers can take 
immediate action to remove the fault and avoid the burden 
of rework and defective items. This objective can only be 
achieved if we have complete, crisp, and authentic informa- 
tion/data about the manufacturing unit under study. Life test- 
ing experiments are very common in the statistical process 
control literature when a random sample of items is selected 
and put on testing for possible failure [2]. The sudden death 
testing 1s repeatedly applied by many parts manufacturers 
with the objective of reduced time of test, see [3]. In this 
technique, a sample of items is first spread in g groups with 
each group having r items. Sudden death testing in groups 
is preferred over the single tester as it reduces the testing 
cost of the product. Sudden death testing technique has 
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been explored by many quality control researchers includ- 
ing [3—6]. It is common to use the Weibull distribution for 
modeling lifetime phenomena. The density function of the 
Weibull distribution is typically written as 


m/(x m1 x \m 
Pa) = = (4) exp|-(=) hoo mand A> 0 (1) 
where À and m are the scale and shape parameters, respec- 
tively. When the shape parameter is one, then this distribu- 
tion tends to form the exponential distribution. The Weibull 
distribution has been explored by many quality control 
researchers including [7] developed life-testing sampling 
plan for the two-parameter Weibull distribution, [8] devel- 
oped bootstrap control chart for Weibull percentiles, [3] 
developed a variables sampling plans for Weibull distrib- 
uted lifetimes under sudden death testing, [9] given progres- 
sively censored reliability sampling plans for the Weibull 
distribution, [10] suggested a group acceptance sampling 
plan for truncated life test having Weibull distribution, [11] 
developed acceptance sampling plans for multi-stage process 
based on a time-truncated test for Weibull distribution. 
Average run length (ARL) is a traditionally used tech- 
nique in the literature of the control charts for evaluating 
the efficiency of the proposed control chart [1]. It is defined 
as the average number of samples when the process is in 
control until the process falsely shows an out-of-control 
signal and is denoted by ARLo, see [12]. To keep the level 
of ARL, at an acceptable level, it should be large enough. 
Another characteristic of ARL is used for the out-of-control 
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process which should be smaller enough for the efficient 
control chart scheme and is denoted by ARL,. ARL for the 
evaluation of a better control chart scheme has been used by 
many authors including [13—24]. 

In practice, however, the observations about the manu- 
facturing process are imprecise, incomplete, uncertain, or 
fuzzy [25]. Fuzzy-based control charts are applied in indus- 
try for the monitoring of the process when some observa- 
tions/parameters are unknown or imprecise. The traditional 
control chart based on the assumption that all observations 
should be determined cannot be applied for the process 
monitoring when some values are fuzzy. Several authors 
contributed in this area by designing a fuzzy control chart 
for various situations, see for example, [26—34, 35]. More 
details can be seen in [36]. 

[37] developed the notion of neutrosophic statistics (NS) 
using the idea of neutrosophic logic was given by Sma- 
randache [38]. The difference between fuzzy logic and 
neutrosophic logic can be seen in [39]. [40-42] discussed 
neutrosophic logic in a variety of fields. The neutrosophic 
statistics is applied when the data in hand are uncertain, 
imprecise, and indeterminate. The neutrosophic statistics is 
more informative than classical statistics. The neutrosophic 
statistics gives the information about the measure of inde- 
terminacy. In addition, it gives the estimated values in the 
indeterminate intervals which are required in uncertain envi- 
ronment, see [37]. Several authors studied the neutrosophic 
statistics including [43] used the neutrosophic interval sta- 
tistical numbers for expressions rock joint roughness coef- 
ficients, [44] studied the anisotropy for neutrosophic num- 
bers of rock joint roughness using the neutrosophic statistics, 
[2] proposed a plan for neutrosophic statistics for testing of 
grouped product using the Weibull distribution, [45] devel- 
oped the monitoring methodology for process variability 
using the neutrosophic statistics, [46] developed a sampling 
plan using the neutrosophic regression estimator. Recently, 
[47] introduced the NS in the area of statistical quality con- 
trol. [45, 48, 49] proposed attribute np, X-bar, and variance 
control charts using NS, respectively. More details on NS- 
based control charts can be seen in [50, 51]. 

[2] designed the sudden death sampling plan using 
the NS. By exploring the literature and to the best of our 
knowledge, there is no work on the control sudden death 
control chart for monitoring the process of having incom- 
plete, imprecise, and vague production data. In this paper, 
we will propose the sudden death control chart for monitor- 
ing the reliability when the lifetime of the product follows 
the Weibull distribution under the uncertainty environment. 
We will present some necessary measures to evaluate the 
performance of the proposed control chart under the neu- 
trosophic statistical interval method (NSIM). We hope that 
the proposed control chart will be more adequate, adequate, 
and informative than the sudden death chart under classical 
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statistics in an uncertain environment. The rest of the paper 
is organized as follows: the design of the proposed chart is 
given in "Design of proposed control chart". The advan- 
tages of the proposed chart are given in "Advantages of the 
proposed control chart". In "Example", an example is given 
and some concluding remarks are given in the last section. 


Design of proposed control chart 


Suppose that Xy =X + ul; I € linf,sup| be a neutrosophic 
random variable having a variable X under classical statis- 
tics (determinate part) and w/;/ € linf,sup| is a indeterminate 
part. Let Xy € La Xy| denote the lifetime of a part, where 
X, is lower value of indeterminacy interval and Xy be the 
upper value of the indeterminacy interval, which follows, the 
neutrosophic Weibull distribution with neutrosophic scale 
parameter Ay € AL. i, | and neutrosophic shape parameter 
my € [m m; |, then the cumulative distribution function can 
be written as: 


F (xy) = 1 — exp(—(Ayay) %); my € [mm], Ay € [A Az] 
(2) 

The operational procedure of the proposed neutrosophic 
control chart for monitoring reliability using sudden death 
testing under Weibull distribution is stated as follows: 

Step-1 Select a random sample of size ny;ny € [n nyl 
having neutrosophic numbers from a manufacturing process 
and distribute ry;řyn € [ri ry| items into gnign € ISL» gul 
groups having ny = ryn- 

Step-2 Perform sudden death testing and detect Y,,,, the 
time to the first failure from the ith group (i = a ee gN) 
and compute teh a Vy = 2o Y NOVN E te vyl and con- 
vert it to VN = VN Ve = [vivi] 

Step-3 Plot vý, on the control chart. Declare the process 


in-control if EOL <7 0c. If 
Vi < LCL sor VN > UCL! 7 then declare the process out- 


of-control. Two control limits, namely LCL? 
eee ee UCL“? 


, 1/3 1/3 1/3 
SUCL,, = UCL, , UCL, , for the proposed control 


chart have been constructed. 

Following Jun et al. [3], the quantity vy € vb. vyl follows 
the neutrosophic Weibull distribution with shape parameter 
2y € ISL» gul and scale parameter ryAy. For developing 
symmetric Shewhart-type control limits, we must transform 
vy € [vi, vyl into a random variable having a neutrosophic 
symmetric distribution. [52] mentioned that transformation 
of vọ = vw * leads to an approximately neutrosophic normal 
distribution with neutrosophic mean 


and 
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ry Ae) PD (an + 1/3 
Hy = (ran) Tulen + 1/3), ime dia (3) 
Ty (gn) 


And neutrosophic variance 


i (raan) Tulen + 2/3) 


2 
vN ` oS - (uyn) > ON E caer | 


(4) 
where I(x) denotes the neutrosophic gamma function. 
Therefore, the control limits for the proposed control chart 

are given as: 


Coe T (en + 1/3) 
Fulgy) 


aa ES 





T- (mn) 
TN (gy) (5) 


— ky 


Py (gy + 1/3) 


LCL? =r a 
N N“ON ly (gn) 





Ty (gy + 2/3 


“Sa ) 


Ty (gy + 2/3) 












—k oe 
E Ty (gn) ( 
(6) 
my \ 1/3 
i= (ryAgn) Tulen + 1/3) 
Ty (gn) 
(rnay Tlen + 2/3) > 
+ ky] —— = lin) 
ly (gn) (7) 
UCL? = nA { Twlgu + 1/3) 
N N*“ON ly (gn) 
P Py (gn +2/3) (a + 2/3) ) 
ly (gn) Ty (gn) 
(8) 


Note here that ky € ky. ky| is neutrosophic control limits 
coefficient and Ap, is the neutrosophic scale parameter of the 
in-control process. 

For the proposed control chart, the process is declared to be 
out-of-control if va (LCLY yor vx JUCL,/ a Therefore, the 


probability that the process is declared as out-of-control when 
the process is actually in-control is given as follows: 


* 1/3 x 1/3 
Posy = P{ ve LCL? doy } + Pf vy > UCLA lon} (9) 


ou 


where 


Pt yi < LCL lon } = PY Ary An vss 2h Ag LCLN |o} 









N 
or 
m Fylen + 1/3 
PL vy < LOL "lon } = oas (2ean B 
N 
Dy (gy + 2/3) (820) 
= ea a a 
Ty (gn) Ty (gy) 
and 


I + 1/3 
Pt vi > UCL)" dow } 4 a aafe a 
N (EN 








Ty (gy + 2/3) Ty (gy + 2/3 


2 
ly (gn) -( Ty (gn) 


So, the average run length (ARL) using the neutrosophic 
statistical interval method (NSIM) is known as the neutro- 
sophic average run length (NARL) is denoted by ARL, which 
was introduced by [45] as follows: 





+ky 


1 
ARLony = ——;ARLon € [ARLo., ARLou] (10) 


outN 
When the process faces a shift, then the NARL is denoted 
by ARL,y, then 
Pan = PY ve (LOLA? | day = coy } + P{vh > UCL? [ayy = cå} 
(11) 


where 


PÍ vi <LCLY don } = P{ 2rydttv( 2ra LCLN’ |A } 








N 
or 
Ty (gyn + 1/3 
P T < LOL" lån } = Goon (2c { T'y (8x + 1/3) 
Dy (gn) 
N Ty (gn) Ty (gy) 

and 
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Ty (en + 1/3) 


P{v* > UCL?|/,} =1-G enad 
{ N i} al N”ON Paw) 








Iy(gy+2/3) [Ty(en +2/3 


Pals) (>52) 


Therefore, the neutrosophic average run length for the 
shifted process is given by 





+ky 


1 


ARL = —$—$ <_< — qm 
P fv (LCL Jain = càon } + P{v* > UCL? |AN = cAgy} 


ARLıy € [ARL;,,ARLiy| 

(12) 
Let doy € ldo dou| denotes the specified NARL. 
The values of ky € ky. kyl will be determined such 
that ARLon € [ARL ARLoy| should be very close to 
don € [dor doy]. The values of ARL,y € [ARL,, ,ARLyy| 
for various values of mp€ [mm] ry € [5,5] 
An € Ags A, | and don € ldo. doy] are determined and are 
shown in Tables 1, 2, 3 and 4. From Tables 1, 2, 3 and 4, 

we note the following trend in neutrosophic parameters. 


1. For all other the same neutrosophic parameters, 
the values of ARL jy € JARL,, ,ARLjy| increase 
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Table 2 The values of NARL when ry € [5,5], my € [0.85, 0.95] and 
various ARL jy, 


ky  [2.470595, 3.771832] [2.507942, 3.831339] [2.520868, 
3.858911] 
Ay __[0.45, 0.55] [0.45, 0.55] [0.45, 0.55] 
gn U1, 3] [1, 3] [1, 3] 
my — [0.85, 0.95] [0.85, 0.95] [0.85, 0.95] 
c ARLy 
1.0 (205.79, 211.75] [315.48, 317.88] [371.48, 387.03] 
1.1 [189.82, 165.16] [290.97, 247.16] [342.61, 300.52] 
1.2 [176.32, 131.88] [270.26, 196.75] [318.22, 238.91] 
1.3  [164.76, 107.4] [252.532 159.75] [297;32; 193:72] 
1.4 [154.73, 88.95] [237.13, 131.9] [279.2, 159.74] 
1.5  [145.94, 74.74] [223.65, 110.5] [263.33, 133.65] 
1.6 [138.18, 63.6] [211.74, 93.74] [249.3, 113.23] 
1.7 = [131.26, 54.71] [201.13, 80.41] [236.8, 97.01] 
1.8  [125.06, 47.54] [191.62, 69.66] [225.6, 83.92] 
1.9 [119.47, 41.66] [183.03, 60.87] [215.49, 73.24] 
2.0 [114.39, 36.8] [175.25, 53.61] [206.31, 64.42] 
Ze [9472272] [145.06, 31.2] [170.76, 37.26] 
3.0 [81.19, 14.38] [124.3, 20.36] [146.31, 24.17] 
4.0  [63.69, 7.80] [97.45, 10.76] [114.68, 12.62] 
5.0  [52.77, 5.05] [80.7, 6.79] [94.96, 7.88] 
6.0 [45.27, 3.66] [69.19, 4.8] [$1.4, 5.51] 


Table 1 The values of NARL when ry €E [5,5], my E [2,2.1] and 
various ARL jy 


ky  [2.746493, 3.548525] [2.848831, 3.636265] [2.925756, 
3.938605] 

Ay  [0.45, 0.55] [0.45, 0.55] [0.45, 0.55] 

en (4, 61 [3, 6] [4, 8] 

my [2, 2.1] [2,2.1] [2, 2.1] 

c ARLy 

1.0  [200.33, 215.34] [306.14, 304.81] [372.56, 374.62] 

1.1 = [181.79, 85.53] [342.51, 118.82] [343.12, 117.25] 

1.2 [106.04, 38.73] [230.25, 52.72] [195.8, 43.99] 

1.3 [62.52, 19.64] [149.42, 26.17] [113.3741931] 

14 [38.84, 11.01] [100.15, 14.33] [69.14, 9.74] 

15 [25.33, 6.74] [69.4, 8.56] [44.22, 5.54] 

16  [17.25, 4.46] [49.55, 5.53] 93,132] 

1.7 = [12.21, 3.16] [36.32, 3.83] [20.44, 2.45] 

1.8  [8.95, 2.39] (27.2152.82] [14.65, 1.85] 

1.9  [6.77, 1.9] [20:92 2.19] [10.83, 1.51] 

2.0 [5.28, 1.59] [16.36, 1.79] [8.23, 1.3] 

2.5 [2.15, 1.06] [6.1, 1.09] [2.96, 1.01] 

3.0 [1.35, 1] [3.09, 1.01] [1.64, 1] 

4.0 [1.02, 1] [1.45, 1] [1.06, 1] 

5.0 [I] [1.09, 1] [1] 

6.0 [1] [1.01, 1] [1] 
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when Ay € AL, A, | increases from Ay €E [0.5, 0.5] to 
An E [1, 1]. 

2. For all other the same neutrosophic parameters, 
the values of ARL,y € |ARL,,,ARLjp| decrease 
when Ay € AL, i, | increases from Ay E [1.5, 1.5] to 
An € [2.0, 2.0]. 

3. From the above observations, we note the increas- 
ing the trend in ARL,y € [ARL ARL] when 
An € |á, A, | <1 and decreasing trend when 
Ay E [Ay AL] > L 

4. From these Tables 1, 2, 3 and 4, it can be seen that the 
value of ARLy decreases as the shifts from 1.0 to 6.0 
increases. 


The following algorithm is used to find the values of 
ky € [kı ky] and ARL,y € [ARL,,,ARLjy|. 

Step-1: determine the suitable range of ky € [kr kul 
and prefix the values of doy € ldo; dou; my € [m,m 
and Ay € |A,, ål. 

Step-2: Determine the values ky € lk. kyl such that 
ARLon = don. During the simulation, we note that sev- 
eral combinations are available which satisfy the given 
condition. Select those values of ky € kL. kul, where 
ARLon € [ARL ARLoy | is very close to doy € dor doy. 
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Table 3 The values of 


NARL when ry € [5, 5], ky [2.366358,3.433075] 

my € [1.45, 1.55] and various 

ARL iy Ax (0.45, 0.55] 
&N [1, 3] 
My [1.45, 1.55] 
c ARLy 
1.0 [207.58, 210.12] 
1.1 [181.24, 140.24] 
1.2 [159.87, 97.47] 
1 [142.42, 70.1] 
1.4 [127.96, 51.91] 
1.5 [115.82, 39.43] 
1.6 [105.52, 30.63] 
t7 [96.68, 24.26] 
1.8 [89.03, 19.56] 
1.9 [82.36, 16.01] 
2.0 [76.49, 13.3] 
25 [55.49, 6.25] 
3.0 [42.71, 3.63] 
4.0 [28.32, 1.85] 
5.0 [20.63, 1.31] 
6.0 [15.96, 1.11] 


Table 4 The values of NARL when ry € [5,5], my € [2.1,2.3] and 


[2.412716, 3.820856] 


[0.45, 0.55] 
[1, 4] 
[1.45, 1.55] 


[304.71, 301.92] 
[266.06, 180.12] 
[234.67, 113.53] 
[209.02, 74.97] 
[187.78, 51.52] 
[169.95, 36.66] 
[154.81, 26.89] 
[141.83, 20.26] 
[130.59, 15.64] 
[120.78, 12.33] 
[112.16, 9.92] 
[81.29, 4.22] 
[62.52, 2.41] 
[41.37, 1.33] 
[30.07, 1.08] 
[23.21, 1.01] 


[2.98 1687, 
3.86765 | 


[0.45, 0.55] 
[2, 4] 
[1.45, 1.55] 


[378.51, 387.3] 
[289.21, 229.92] 
[226.43, 144.18] 
[180.95, 94.7] 
[147.15, 64.72] 
[121.49, 45.79] 
[101.63, 33.39] 
[86.01, 25.01] 
[73.54, 19.19] 
[63.46, 15.04] 
[55.21, 12.02] 
[30.41, 4.96] 
[18.94, 2.73] 
[9.29, 1.42] 
[5.56, 1.11] 
[3.78, 1.02] 


various ARL jy, 


kx  [2.785936, 3.654638] [2.39096, 3.434913]  [2.418024, 


3.773643] 
An  [0.55, 0.65] [0.55, 0.65] [0.55, 0.65] 
gn [2,4] [1, 3] [1, 4] 
my [2.1, 2.3] (2:1; 2.3] [2.1, 2.3] 
c ARLy 
1.0  [205.3, 206.33] [300.76, 302.71] [379.18, 372.29] 
1.1  [139.68, 97.49] [251.72, 165.29] (317.96, 172.57] 
1.2 [98.53, 50.47] [210.24, 96.18] [265.58, 87.46] 
1.3 [71.68, 28.28] [177.83, 59.11] [224.61, 47.88] 
14 [53.54, 16.99] [152.27, 38.09] [192.31, 28.04] 
1.5  [40.93, 10.86] [131.8, 25.59] [166.44, 17.44] 
16 [31.93, 7.33] [115.16, 17.85] [145.41, 11.44] 
1.7 = [25.36, 5.21] [101.45, 12.87] [128.09, 7.88] 
1.8  [20.46, 3.87] [90.03, 9.57] [113.66, 5.68] 
1.9 [16.76, 3] [80.42, 7.31] [101.51, 4.25] 
2.0  [13.9, 2.41] [72.26, 5.73] [91.2, 3.31] 
2.5 [6.43, 1.26] [45.42, 2.33] [57.27, 1.47] 
3.0  [3.65, 1.03] [31.13, 1.42] [39.21, 1.09] 
4.0 [1.79, 1] [17.24, 1.03] [21.66, 1] 
530 [125.1] [10.98, 1] [13.75, 1] 
6.0 [1.07, 1] [7.66, 1] [9.54, 1] 


Step-3: Using the selected values of ky € kL kyl, deter- 
mine the values of ARL,y € [ARL L, ARL;y] for various c. 


Advantages of the proposed control chart 


In this section, we will discuss the advantages of the pro- 
posed control chart under neutrosophic statistics over the 
existing control chart under classical statistics proposed by 
[53]. [43, 44] mentioned that a method which provides the 
output in indeterminate interval under uncertainty is known 
as the most an effective and adequate method. To compare 
the proposed chart with [53], we will fix the same values of 
the control chart parameters. The values of NARL of the 
proposed control chart and [53] chart are shown in Table 5 
for various values of control chart parameters. 

From Table 5, it can be noted that all values of [53] 
chart are within the indeterminacy interval of the pro- 
posed control chart. The NARL values of the proposed 
control chart become the same as in [53] if no uncertain 
observation is recorded in the data. For example, when 
c=1.1 and Ay € [2,2], the indeterminacy interval of 
NARL for the proposed chart is from 172nd sample to 
the 317th sample. It means that one can expect that the 
proposed control chart indicates that on average, the pro- 
cess will be shifted between 172nd to 317th samples. On 
the other hand, [53] indicates that the first out-of-control 
will be detected at the 317th sample. From this study, it 
can be seen that the proposed chart gives smaller values 
of NARL as compared to the existing chart. In addition, 
from this comparison, it is clear that the proposed control 
chart gives the results in indeterminate interval; therefore, 
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Table 5 The comparison in NARL 
c  NARL 
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ARL 


Proposed chart when J = 1.0 Proposed chart when A= 1.5 Proposed chart when 4 = 2.0 Existing chart Existing chart Existing chart 


whenA=1.0 whenA=1.5 when A= 2.0 


1.0 
1.1 
1.2 
1.3 
1.4 
1.5 


[371.48, 387.03] 
[342.61, 300.52] 
[318.22, 238.91] 
[297.32, 193.72] 


[279.2, 159.74] 


[263.33, 133.65] 


[378.51, 387.3] 


[289.21, 229.92] 
[226.43, 144.18] 


[180.95, 94.7] 
[147.15, 64.72] 
[121.49, 45.79] 


1.6 [249.3, 113.23] [101.63, 33.39] 
1.7 [236.8, 97.01] [86.01, 25.01] 
1.8 [225.6, 83.92] [73.54, 19.19] 
1.9 [215.49, 73.24] [63.46, 15.04] 
2.0 [206.31, 64.42] [55.21, 12.02] 
2.5 [170.76, 37.26] [30.41, 4.96] 
3.0 [146.31, 24.17] [18.94,2.73] 
4.0 [114.68, 12.62] [9.29, 1.42] 
5.0 [94.96, 7.88] [5.56, 1.11] 
6.0 [81.4, 5.51] [3.78, 1.02] 


itis flexible, informative, and adequate and concurs with 
the theory of [43, 44]. Hence, the proposed control chart 
under neutrosophic statistics is a better alternative of [53] 
chart under classical statistics when there is Neutrosophy 
in the data. 

Now, we present the efficiency of the proposed control 
chart over [53] using the simulated data. Using the neu- 
trosophic Weibull distribution, the neutrosophic data were 
generated in which the first 20 observations belong to the 
in-control process while the next 20 observations from the 
shifted process with c=2.0. Now according to Table 3, 
ARL,y € [55.21, 12.02] when ky €E [2.981687, 3.86765] 
, An E [0.45, 0.55], en € [2,4] and my €[1.45,1.55]. 
According to the proposed chart, it can be expected a 
shift between 55th sample and 12th sample. The values 
of statistic vý € lv, v | are plotted in Fig. | for the pro- 
posed chart and in Fig. 2 for [53] chart under classical 
Statistics. We note from Fig. 1 that the 29th sample is 
at the neutrosophic lower limit which indicates that the 
process has shifted. From the figure, it can be seen that 
all plotting values are within the control limits and do 
not show any shift in the process. From Figs. 1, 2, it is 
clear that the proposed control chart detects a shift in the 
process at the 29th sample while the existing control chart 
indicates that the process is an in-control state. Figure 2 
also shows several values of vý € lve, | within the inde- 
terminacy interval. From this comparison, it is quite clear 
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[379.18, 372.29] 371.48 378.51 379.18 
[317.96, 172.57] 342.61 289.21 317.96 
[265.58, 87.46] 318.22 226.43 265.58 
[224.61, 47.88] 297.32 180.95 224.61 
[192.31, 28.04] 279.2 147.15 192.31 
[166.44, 17.44] 263.33 121.49 166.44 
[145.41, 11.44] 249.3 101.63 145.41 
[128.09, 7.88] 236.8 86.01 128.09 
[113.66, 5.68] 225.6 73.54 113.66 
[101.51, 4.25] 215.49 63.46 101.51 
(91.2, 3.31] 206.31 55.21 91.2 
[57.27, 1.47] 170.76 30.41 3/:42/ 
[39.21, 1.09] 146.31 18.94 39.21 
[21.66, 1] 114.68 9.29 21.66 
[13.75, 1] 94.96 5.56 13.75 
[9.54, 1] 81.4 3.78 9.54 
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Fig. 1 The proposed chart for the simulated data 


that the proposed chart has the ability to detect a shift in 
the process. 


Example 


Suppose that a ball bearing company is interested to apply 
the proposed control chart for the monitoring of the ball 
bearing reliability. The company found that the reliabil- 
ity data follow the neutrosophic Weibull distribution with 
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Fig. 2 [53] chart for the simulated data 


mye[2.1,2.3] and Aye[0.55,0.65]. For this experiment, let 
ryEel5,5] and ARL e[370, 370]. For these parameters, 
enel4, 4] from Table 4. The experimenter will select a 
random sample of 20 balls and distribute them into 4-sub- 
group size. During the experiment, it is found that the 
reliability values of ball bearing are neutrosophic numbers 
rather than the exact numbers. Therefore, the monitoring 
of the process will be done using the proposed control 
chart. The neutrosophic data are shown in Table 6. 


Table6 The my = [2.1, 2.3], 


The values of statistic viel|v*, v | are plotted in Fig. 3 
for the proposed control chart and in Fig. 4 for the exist- 
ing control chart proposed by [53] under classical statistics. 
From Figs. 3, 4, we note that the proposed control chart 
provides the values in indeterminacy interval rather than the 
exact values. In addition, the proposed control chart shows 
that three plotting statistics are near the control limit which 
should be a problem in the process. 


Conclusion 


In this article, a neutrosophic control chart for monitoring 
the reliability using sudden death testing under Weibull 
distribution has been presented. The NSIM has been used 
for the required measures to apply the proposed control 
chart. A simulation study was made to show the compara- 
tive efficiency of the proposed chart for the uncertainty 
environment. From the comparative study, it was noted 
that the proposed chart is efficient than the existing chart. 
In addition, on comparing the proposed chart, it can be 
observed that it is adequate, more flexible, and more effi- 
cient for use in an uncertain environment. In nutshell, the 
proposed chart is a valuable addition in the toolkit of the 
quality control personnel when there are unclear, uncer- 
tain, vague, or fuzzy observations in the sample. The pro- 
posed chart can be extended for the cost models or some 
other sampling schemes as future research. 


V*[N] 


Sr.# Sample 

An = [0.55, 0.65] in-control 

process data generated from 1 [0.0178, 0.0178] 

Weibull distribution 2 [0.4935, 0.4935] 
3 [0.1535, 0.1535] 
4 [0.4487, 0.4636] 
5 [0.5741, 0.5741] 
6 [0.3445, 0.3445] 
7 [0.5875, 0.5875] 
8 [0.3616, 0.3616] 
9 [0.4093, 0.3518] 
10 [0.1626, 0.1626] 
11 [0.3981, 0.3981] 
12 [0.5814, 0.5814] 
13 [0.1800, 0.1800] 
14 [0.5739, 0.5739] 
15 [0.2634, 0.2634] 
16 [0.3777, 0.3777] 
17 [0.3264, 0.1707] 
18 [0.4717, 0.4717] 
19 [0.5016, 0.7292] 
20 [0.6478, 0.6478] 


[0.5117, 0.6457] 
[0.6793, 0.6793] 
[0.4894, 0.4894] 
[0.7130, 0.7130] 
[0.7975, 0.7975] 
[0.5629, 0.5629] 
[0.8473, 0.5624] 
[0.4614, 0.7704] 
[0.7233, 0.7179] 
[0.4519, 0.4519] 
[0.5113, 0.5113] 
[0.6823, 0.6823] 
[0.3272, 0.3272] 
[0.6138, 0.6138] 
0.4254, 0.4254] 
[0.4797, 0.4797] 
[0.4104, 0.4104] 
[0.8010, 0.8010] 
[0.5697, 0.5697] 
[0.6580, 0.6580] 


[0.7653, 0.7653] 
[0.9169, 0.9169] 
[0.6024, 0.6024] 
[0.7811, 0.7811] 
[0.8576, 0.8576] 
[0.5748, 0.5748] 
[1.0369, 1.0369] 
[0.6973, 0.8543] 
[0.9794, 0.9794] 
[0.7953, 0.7953] 
[0.8941, 0.8941] 
[0.8180, 0.8180] 
[0.5806, 0.5806] 
[0.6609, 0.6609] 
[0.5915, 0.5915] 
[0.6254, 0.6254] 
[0.6024, 0.6024] 
[0.9415, 1.2846] 
[0.6646, 0.6646] 
[0.8051, 0.8051] 


[0.8343, 0.8734] 
[1.0457, 1.0457] 
[0.6175, 0.6175] 
[0.8199, 0.8199] 
[0.8578, 0.8578] 
[0.6710, 0.8392] 
[1.2262, 1.2262] 
[0.9236, 0.9236] 
[1.0072, 1.0072] 
[0.8459, 0.8459] 
[0.9158, 0.9158] 
[0.8316, 0.9851] 
[0.9519, 0.9519] 
[0.8334, 0.8334] 
[0.8646, 0.8646] 
[0.7428, 0.9642] 
[0.6308, 0.6308] 
[1.0895, 1.0895] 
[0.7038, 0.7038] 
[1.0139, 1.0139] 


[1.1012, 1.1599] 
[1.3526, 1.3526] 
[0.9119, 0.9119] 
[1.1963, 1.1989] 
[1.2967, 1.2967] 
[0.9750, 1.0636] 
[1.54460, 1.483] 
[1.1316, 1.2794] 
[1.3607, 1.3531] 
[1.1107, 1.1107] 
[1.2266, 1.2266] 
[1.2338, 1.3027] 
[1.0673, 1.0673] 
[1.1533, 1.1533] 
[1.0378, 1.0378] 
[1.0106, 1.1369] 
[0.9151, 0.8953] 
[1.4138, 1.5654] 
[1.0623, 1.1382] 
[1.3211, 1.3211] 
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Fig.3 The proposed chart for ball bearing data 
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Fig.4 The existing chart for the ball bearing data 
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